Statistical theory of monodisperse linear polymers of finite semifiexibility: 

Bethe-Peierls approximation 
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We employ the Bethe-Peierls technique to obtain the self-consistent thermodynamic description 
for the system of linear semiflexible polymers of equal length. To our best knowledge, this is 
done rigorously for the first time. We show that the theory is in quantitative agreement with 
the Sanchez-Lacombe theory. The overall configurational entropy expression contains the athermal 
and the semifiexibility contributions and the correction due to the attractive interactions between 
monomers. Based on the paper by Wolfgardt et al. [Phys. Rev. E 54, 1535 (1996)], we conclude 
that the present theory captures the first two contributions better than the Gibbs-DiMarzio theory. 
The entropies of the states corresponding to high densities, obtained at low temperatures, high 
pressures and high molecular weights, extrapolate to negative values. This feature, known also as 
the Kauzmann paradox, is the thermodynamic manifestation of the glass transition. An investigation 
of the configurational specific heat dependence on molecular weight is also provided. 



I. INTRODUCTION 

In 1948, Kauzmanrti recognized a peculiar feature of 
the thermodynamic behavior of liquids in the vicinity of 
glass transition. He found that the configurational en- 
tropy drops very rapidly with temperature decrease and 
does not seem to vanish close to absolute zero. It looks 
like the configurational entropy would instead drop to 
zero at a, positive temperature, 7k, and, hence, something 
should happen to the system in order to avoid the en- 
tropy becoming negative at temperatures below Ik- The 
problem has attracted much theoretical attention after 
the work by Gibbs and Di Marzic 2 . They demonstrated 
by means of a statistical calculation that the configura- 
tional entropy of a disordered polymer liquid extrapolates 
to the negative values at low temperatures. This vio- 
lates the Nernst postulate and makes the problem very 
puzzling. To resolve this problem, they proposed that a 
second-order transition to a unique state, the ideal glass, 
would occur at T = Tk- Gujrati and Goldstein 3 crit- 
icized the Gibbs-DiMarzio theory for the absence of a 
crystalline state. But, recently, Corsi and Gujrati^ pre- 
sented a calculation where, besides the disordered liquid 
phase, the crystalline phase was captured and where the 
configurational entropy of the liquid phase would become 
negative at low temperatures as well. It was pointed out 
that, provided the stable crystalline state exists below the 
melting temperature, Tm > Tk, the unusual ideal-glass 
transition in the metastable liquid region by no means 
should violate the thermodynamic laws, see also^. 

The issue of the crystalline state existence for complex 
polymer topologies is a very challenging theoretical prob- 
lem. In the course of polymerization chains of progres- 
sively larger sizes are formed, which causes the increase 
in viscosity and eventually leads to the vitrification of the 
system. In the case when the polymerization coordina- 
tion is larger than two, the polymers can be formed with 
branches extending in an irregular manner. This makes 
it difficult if not impossible to construct an ordered state 
in the system. 



In the present paper we study a model of compress- 
ible polymer melt by means of a new treatment based 
on the Bethe-Peierls approximation. As in the Gibbs- 
DiMarzio theory, our calculation does not capture the 
ground state. The latter fact is, of course, a disadvan- 
tage of the method, but the level of prediction is quali- 
tatively the same as in the more accurate theory^. We 
show that the present simple theory enjoys a quantitative 
agreement with experiment and can be used in connec- 
tion with the work by Wolfgardt et al&. 



We consider the semiflexible linear x-mers distributed 
on a lattice. The penalty for a chain bend is intro- 
duced by means of the Flory's flexibility 7 , /. We also 
take into account the interactions among monomers as- 
signing the energy, V < (attractive interaction), to 
the nearest-neighbor pairs of occupied sites. The Bethe- 
Peierls approximation^ was originally developed for sin- 
gle site species. An extension of the method to diatomic 
species was provided by Chang^. The aim of the present 
paper is to extend the Chang's scheme further to poly- 
atomic species. A geometrical interpretation of polymer 
semifiexibility, consistent with the Flory's definition of/, 
can also be provided, which is demonstrated in the fol- 
lowing. The entropy expression is found to reduce to the 
Flory-Huggins formic in the appropriate limit: V = 
and z — * oo, where z is the lattice coordination number. 



A significant improvement to the Flory-Huggins the- 
ory is the lattice cluster theory (LCT) developed by 
Freed and co-workers (seeAi and references there in) . The 
LCT not only provided a rigorous and systematic clus- 
ter expansion method for generating corrections to the 
Flory-Huggins theory, but also predicted the influence 
of monomer structure on thermodynamic properties an- 
alytically. However, the formalism is rather complicated, 
which impedes its wide use. In analogy to the Gibbs- 
DiMarzio theory, the LCT configurational entropy also 
extrapolates to zero at a positive temperature^ 2 -. 
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II. BETHE-PEIERLS APPROXIMATION 

The main idea of the method is to calculate the grand 
partition function (GPF) for configurations of an ag- 
gregate composed of a molecule occupying a lattice site 
(internal site) and its z neighbors (surface sites) in the 
framework of a mean-field theory: An effective field is 
superimposed on the surface molecules by molecules out- 
side the aggregate. It is determined making use of the 
condition of equal probability for the internal molecule 
as well as for any molecule from the aggregate surface 
to be in a particular configuration. After determination 
of the GPF other quantities such as energy (enthalpy at 
constant pressure), specific heat, etc. follow at once. 

The effective interaction parameter, V, is defined as 
the energy of interaction between two occupied sites 
that are nearest-neighbors to each other irrespective of 
whether they are connected by a bond or not. To ac- 
count for the presence of free volume, the lattice-hole 
model is employed in which monomers and holes (empty 
lattice sites) populate a lattice with the occupancy con- 
trolled by the system chemical potential, fi. The polymer 
scmiflexibility is incorporated through U / being the ener- 
getic penalty for creation of a chain bend. Based on this 
information, a general form of the GPF can be given: 

Z = J2n(n,X,N g ,N) V x \ n w^, (1) 

where r) = exp(— V/k^T) is the Boltzmann weight for the 
interactions between occupied sites, A = exp(/i//cBT) is 
the absolute activity, Wf — exp(— Vf/k^T) is the Boltz- 
mann weight for the flex energy, n is the total number 
of occupied sites, X is the number of nearest-neighbor 
pairs of occupied sites, N g is the total number of the 
gauche configurations, Q is the number of arrangements 
on a lattice of N sites that have the same n, X, and N g 
numbers. 

For comparison purposes, it is instructive to begin by 
presenting the original Bethe-Peierls technique. Next, we 
are going to present a modification due to the chain con- 
nectivity. At first, the aggregate is pictured as a site sur- 
rounded by z surface sites. A set of numbers is defined: 
{0i} = 0q, 6\, 9 Z , which are or 1 according as the cor- 
responding sites are empty or occupied, where the index 
i = is reserved for the internal site. A configuration 
with the internal site being unoccupied (0; 9\, ...,9 Z ) con- 
tributes to the GPF as the sum over all arrangements on 
the surface sites J2{8t } \ 6l+ '" +e * ip(0i, 6 Z ), where tp is 
the field that couples the aggregate configurations to the 
uniform external distribution. Presuming there are no 
sites occupied by the segments of the same chain among 
the surface sites, two types of correlations generated by a 
local configuration can be identified: the occupancy of a 
surface site 1) means the possibility of the bond presence 
between the given site and an external nearest-neighbor 
site, 2) induces the mer-mer interaction with occupied ex- 
ternal sites. An empty surface site does not induce any 
correlations to the exterior. These facts can be generally 



taken into account assuming that ip is given by 
^ = c £0i+...+^ 

where c is some constant. From this follows that for the 
configuration (0;#i, ...,9 Z ) the GPF is given by 

Y,c(\0 9l+ - +e > =c(l + XO z - (2) 

To study the linear chains, we consider the aggregate 
composed of x sites and their z' — 2(z — 1) + (x — 2)(z — 2) 
next neighbors. The modification was inspired by the 
papers published quite some time ago2J^, which makes 
us a little unsure about the novelty of the result, but we 
are unaware of any related work. 

The GPF for the configuration (1, 1; 9 x +i, 9 x +z>) 
can be represented in the following general form: 

£ Y. c.\'rTH>^F B+ "' +Trt , (3) 

where the first sum is over configurations permitted by 
the linear chain topology, the second sum is over config- 
urations on the surface sites, the factor i] x ~ 1 takes into 
account the interactions among the segments of the chain 
inside the aggregate, and J2 3 9 = 9x + ■■+9j-i+9j+x + 9 x . 
In order to calculate ([3]) , we propose a recursive method 
analogous to the Bethe lattice technique^. Accordingly, 
the aggregate composed of a site and its z neighboring 
sites corresponds to the center of the Bethe lattice, an in- 
finite Cayley tree, which is a connected dendritic struc- 
ture with constant coordination, z, and no loops. As 
before, the configurations with the central site being un- 
occupied are included in the GPF by the term c(l + A£) z , 
where £ represents the effective field for the case when no 
bond is present inside the aggregate. For the case when 
the central site is occupied by an end-point or by a middle 
group, we define Xk to account for correlations induced by 
the chain connectivity; where k is the number of bonds 
in the chain piece extending beyond the surface site of 
the aggregate. The GPF can be formulated as follows: 

Z = c(l + Xif + cz(l + AT?0" -1 A a »7Xx-i 

x-2 

+c|r/(l + A^r^Ay^TxfcX.-fc^, (4) 

k=l 

where the factor z in the second term is the number of 
orientations for the bond extending from the end-point 
occupying the internal site of the aggregate; in the third 
term, the factor 

r f =w f (z-2) + l (5) 

takes into account the semiflcxibility. The latter is de- 
rived making use of the assumption that the number of 
trans configurations in the aggregate is equal to z/2. In 
the case when z is an even number and the surface sites, 
situated around the internal (central) site, are equally 
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spaced among each other, z/2 is easily identified as the 
number configurations of two bonds with the 180 degree 
angle between them. Note that this naive picture cannot 
be used for the case when z is an odd number and the 
surface sites are equally spaced around the internal site, 
since a construction of trans configurations with the 180 
degree angle between the bonds is not possible in 2-D 
space (we are not sure about the 3-D case). Proceeding 
further with our description of we justify the statis- 
tical weight zrf/2 by means of the relation 



z(z-l) 



where z{z — l)/2 is the total number of possible config- 
urations for two bonds inside the aggregate. In order to 
evaluate the GPF, we obtain the Behte lattice recursive 
relations for Xk- 

xi = c'ii + x^y- 1 , 

X 2 = c'r f \r)(l + \r)(,) z - 2 xu 



Xk = c'r f Xr)(l + \r)O z Xk-i, 

where the factor rj and the power z — 2 are due to the 
fact that there are z — 2 local gauche conformations out 
of z — 1 possible for a middle segment located anywhere 
except the center of the Bethe lattice, c' is a scaling factor 
measuring the distance from the Bethe lattice center to 
be evaluated in the following. The summation in the 
third term of (j3J) can now be carried out to yield 

x—l 

E»Xx- t - 2 = (x - 2)(c'r / Ar 7 )*- 3 (l + \ v tf*- 3 X'-Q '*? 

k=2 

Finally, we find 

Z = c(l + \§ z + af x - 1 ^-r B f - 2 xX B rf- 1 (l + A»?0*'- ( 6 ) 
2 3 

Therefore, the function ip for this configuration ap- 
pears to be 



where 7 is the number of ways per site the polymer chain 
can be arranged on otherwise empty lattice. We find 
that 7flox = §(z — l) 21-1 for completely flexible chains, 
7rigid = § for rigid rods, and 7 flox = §[w f (z - 2) + lf^ 1 
for semiflexible chains. The GPF for this configuration 
will be 

c xl x\ x i 1 x - 1 {l + x^y' , 

where the constants c x and £ x can be related to c and 
£, respectively: We write the GPF for the aggregate 
with the empty internal site and one empty surface 
site, and compare it with the GPF written for the ag- 
gregate with two empty sites being internal. The for- 
mer is obtained summing over all d's on z — 1 sur- 
face sites, which results J2{9 Z } C ( X O 02+ ' " +9z ■ The lat " 
ter is evaluated by summing over 2{z — 1) surface sites: 

'E{8 t ,e>} c 2( X &) e2+ '" +e * +e ' 2+ '" +e ' z - Both partition func- 
tions should be the same, hence the relations £ = £2 and 
c = 02(1 + A£) 2_1 are valid, as found by performing the 
summation over various configurations 9' 2 ,...,0' z . This 
idea can be applied for aggregates of larger sizes. Our 
suggestion: since the chains consist of x — 1 bonds each, 
the relation should be 



£ = £ x and 



c*(l + A£) 



(z-l)(x-l) 



(7) 



c x ix\ x V x - l (\r]Z x y*+i 



+•••+«„ 



which yields 

z = C (i + ao z + c^Av^a + ^iy\ (8) 

which is identical to ((6]) provided c' — (1 + A£)~^ _1 ). 

Having established the GPF, we now turn to the de- 
tailed study of the equal probability condition for absorp- 
tion at the internal site and surface sites. Note also that 
this probability is the fractional coverage of the lattice 
9 = n/N . The following pair of equations 

» /. Y. E <>AV ux,,i,^ " 

= z-^xjxX^-^i + x-nO*' (9) 

and 



(zz)- i cj2(&i + - + o z m) ei+ - +e ' 

+{zz x )- l c x J2 zZ (i_L^_Li + H 6 + - + Yj X e ^ xx ^ 1 ^ Xl ^ T ' 2e+ '" + ^ e 

{zZ)- l ^{l + AO 2 + {zZx)- l c xl X x rf- 1 U{x - 1)(1 + Xrj£) z ' + £^(1 + XrfiA , 



(10) 



are obtained in analogy with the corresponding equations from Ref. 9. Combining Eqs. (19)) and (|10p . one finds the 
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equations of state: 



zx e{\ + rje) 
V l + e 
xe (1 + e) x - 2+z ' 
z' (1 + -qeY'- 1 ' 



(11) 
(12) 



where e — A£. The solution of the quadratic equation 
(fTTj) is given by 



e(9) =e+ = 



6z' - (l-9)zx + D 
2zxr\{\ - 9) : 



(13) 



where 



D = y/\6z> - (1 - 8)zx] 2 + Az'zxri6{\ - 9). 

Another solution with the negative sign before D in [)13p 
is physically irrelevant, since £_ = —1 in the athermal 
limit rj — 1, while e + = z'0/zx(l — 9) > in this case. 

There is another method to derive (fTT|) and (fT2")k The 
occupation is obtained as the ratio of the term corre- 
sponding to the occupied internal site to the total parti- 
tion function: 



8 = Z- l c-r x r 2 x\ x ^- x {\ + AO~ (a!_1)( * _1) (l + XvO Z ' 
2 J 

(14) 

This allows us to derive (|12[) : 



- x - 2 xX x ri x - 1 = 



2 f 



(1 + K) 



x—l+z' 



1 



d + x v o z ' ' 



(15) 



It should also be considered that the probability to find 
the internal site of the aggregate unoccupied can be cal- 
culated by two methods: 



and 



Z- l c{(\ + AC) 2 " 1 +(z- l)AVl + At?0*" 2 X*-i 



z-2 



r/(l + AT ? 0*- 3 AVX)xfcXx-*-i} > 



fc=2 



where the quantity in the curly brackets is the partial 
partition function for the surface z sites calculated under 
the condition that the internal site is empty. Thus, with 
aid of (fT5|) . we obtain 

z' 9 (1 + A0 Z 



(i + Ae) z = (i + Ao z - 1 + 



zxl-e(l + AryO ' 



which can be reduced to (|11[) . 

To this end we remark that one can also derive 
these equations solely from the kinetic considerations 
for molecular absorption/desorption process^. Here we 
provide only the derivation of the first equation to save 
space. In the aggregate represented by a site surrounded 
by z sites, the probability for a configuration (0; 9\, ...,9 Z ) 



to occur is assumed to be proportional to 6 ' 1+ - +9z , 
where 4> is the quantity to be determined further. Then, 
the probability for the internal site to be unoccupied 
should be proportional to Yl{e i } ( t> ei+ '" +6z = (1 + 4>Y ■ 
Hence, the probability for one of its z neighbors to be 
occupied is given by 



(1 



\z-l 



(1 



(1 + 4>Y 



1 



(16) 



For the aggregate enclosing a chain with its neighbor- 
hood, the probability for any of the surface sites to be 
occupied, provided all internal sites are empty, is calcu- 
lated in the analogous manner and is given by 



(17) 



1 



where <f> is the quantity similar in nature to <fi. Comparing 
(fl"6|) and (fT7"|) . one gets <f> — (j). The probabilities for 
various states of occupation of the surface sites, provided 
that all x internal sites are occupied, are proportional to 
(jl£yx+i+---+6 x+S! > ^ wn i cn gives 



1+T]<t> 



for the probability of a surface site to be occupied when 
the internal sites are also occupied. The condition is 



= (1- 



'l + (t> 
2{x - 1) 



zx 1 + r](f> 



where [2{x — l)/zx]9 is the ratio of number of bonds 
2(x — 1)9 /zx per lattice site to the total number of lattice 
bonds z/2 per lattice site representing the probability 
to have a bond extending from an occupied site. The 
equation pip follows immediately taking 



III. THERMODYNAMICS 

Further analysis of |T]) can be carried out to derive the 
configurational entropy expression. First, we write the 
partition function as 



Z = J2 nXnr l X Wf 



_ e N(0g(9)+6 In \+6 x In r)+fl g ln Wj ) 

where 9x — X/N, and 9 g = N g /N . The minimal free 
energy can by achieved applying the condition: 



d_ 

09 



[9g{9) +9\n\ + 9 x In i] + 9 g In w f ] = 0. 
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Thus, the configurational entropy is given by 

S N /k B = [0g(0)Y 

= - J In Xd9 - Ox In J] - 6 g In w / , (18) 

where Sjy is the entropy per lattice site. Note, g*(0) = 0, 
9x(0) = 0, and 9 (O) = 0. It is appropriate to define the 
entropy per monomer: 

S n = Sn/0. 

After integration of (fT5)l . the entropy expression adopts 
the form 

S n /k B = -\n( 1 x7 1 x - 1 )--\n6- < ^—^\n{l-e) 

X XV 



(x-i)(i-e) 

xO 

x x 



ln[(l + £ )(l-0)] 

ln(l - 9) (19) 



-E n /k B T, 



x20 ln \l 



where E n /k B T = — (9xlnrj + 9 g \nwf)/9 is the system 
internal energy per monomer. 

The contact density 9x is evaluated making use of the 
lattice topological relations: 

z9 = 29 X +9 X ,, 
z{l-6) = 26x«+6x', 

where 0x> is the density of pairs formed by an empty site 
adjacent to an occupied and 9x" is the density of nearest- 
neighbor pairs formed by empty sites. We calculate 9x" 
from the fact that the probability of having an empty site 
being a neighbor to another empty site is, in accordance 
with (fT!)]) , equal to 1/(1 + e), and taking into account 
that there are on average z/2 lattice bonds per site: 



7X" = ~ 



* (l-g) 

2 1 + e ' 



Then, we have 



'i x = zb 



z z (1 



2 2 1 + e 
The density of the gauche bonds 9 g is given by 

(z - 2)w f 



-(x-2) 



l + (z- 2) WJ ' 



(20) 



(21) 



where 9/x is the number of polymers, {x — 2) is the num- 
ber of middle groups in each polymer, and the factor 



f=(z-2)w f /[l + (z-2)w f ] 



(22) 



is the probability to form the gauche configuration de- 
rived by Floryi. The latter can be obtained from our 



theory as well: The density of gauche bonds 9 g is the 
ratio of 

z(z - 1) z _ z(z - 2) 
2 2 ~ 2 

terms corresponding to the gauche configurations in (|4|) 
to the total the GPF: 

_ (z/2)(z - 2)w f (x - 2)X x i] x ' 1 (l + AfJ-C*- 1 )^- 1 )^ + Xr^) z ' 



which, by use of Eq. (|14p . can be written in the more 
compact form: 

X Tf 

This completes our proof of (f2"Tj) . 

The equilibrium properties of the system can be de- 
duced from the thermodynamic identity 



fi — H n — TS„ 



(23) 



where the configurational enthalpy per monomer, H n = 
E n +Pvo 1 9, with vq being the volume of a lattice site. By 
aid of Eq. (fT2"| the density dependence of the chemical 
potential, [i = k B Thi\, arises straightforwardly: 



M = 



^^-\n(~/xrj x 1 ) + -^-^ln 



+ *5l( x _l)l n (l +e )+/l n f_L±! 
x \ 1 + rje 



Performing further algebraic manipulations, we derive 
Pv /k B T= f]n(i + e ) + iZL?ln(i-e). (24) 



The equation (|24[) serves as the implicit equation for de- 
termination of 9. 

Making use of the relations for the densities, Eq. |T 
can be written in a slightly different form 

S n /k B = -ln- + ln(l - 9) 

x 2 x x a 



x6 

x x 

i. 1 

X \ z' J X 
Ix - 9 h )\vLr\ 



dn[(l+e)(l-i 



1+e 



(25) 



z . / zxe\ z l , /l 

— In M In - 

x \ z' J x29 VI 



x-2 



[/m/ + (l-/)ln(l-/) 



-/ln(z-2)], 

where 9b is the bond density given by 9{, = (x — 1)9 /x. 

The equations (|24|1 and (fT3|) are to be solved numer- 
ically for 9 as a function of T at constant P. We then 



6 



0.98 















* 




0.96 




















"^j. , P* = 0.25 










■ 


LDPE 






▲ 


PnBMA 






* 


PIB 




0.92 


O 


PS 






A 


PoMS 


a X 




□ 


PcHMA 






* 


PMMA 




0.90 




PDMS 










A 



0.5 0.6 0.7 0.8 

T* 



FIG. 1: Test of Eq. (|24[) against experimental data for various 
polymer liquids. 

substitute 9 and e(6) calculated from (fT"3]) into (|25p to 
obtain the entropy at a constant pressure. Finally, con- 
figurational heat capacity per occupied site, C p , can be 
computed from 

C p = {dH n /dT) p 

or equally from 

C v = T(dS n /dT) p . 

According to our numerical computation, both identities 
produce the same specific heat, which indicates that the 
present theory is thermodynamically self-consistent. 

IV. COMPARISON WITH OTHER RESULTS 

In the present section, we analyze the main charac- 
teristics of the thermodynamic functions (|24| and (125]) 
in comparison with data available from experiments and 
simulations. 

Presented as early as in 1976, the lattice fluid (LF) 
theory developed by Sanchez and Lacombei£, due to its 
simplicity and high level of prediction, became a founda- 
tion of many successional theoretical and experimental 
investigations of polymer systems 17 . In analogy with the 
LF theory, the equation of state (f24|) can be written in 
terms of the reduced temperature T* — —2ksT/zV and 
pressure P* = —2vqT/zV. Figure [T] shows the exper- 
imental data for density, 6, vs. reduced temperature, 
T* , obtained from Fig. 3 of Ref. 18 for various liquids. 
The line is the theoretical isobar (T* = 0.25) computed 
applying the limit x — + oo for z = 20. Based on this 
we conclude that the present theory adequately describes 
the relationship of density, temperature, and pressure for 
pure polymer fluid. 



The important difference between the theory by 
Sanchez and Lacombe and the present theory is in en- 
tropy. The LF entropy is independent of temperature, 
while (|2"5|) exhibits a temperature dependence. This is a 
consequence of consideration of the configurational dis- 
tribution of the polymer system influenced by mer-mer 
interaction energy and polymer semiflexibility, which is 
absent in the LF theory. 

At this point we would like to note that the results from 
the Monte Carlo (MC) simulations by Wolfgart et al~ 
could be used in support of the present theory. These MC 
simulations revealed that, for the case of semiflcxible non- 
interacting polymer, the Milchev's entropy expression 19 

^ /, 1 , z his hi# (1 — 9) , / ox ,„„x 

S M k B = -ln- + ^—^ In 1-0 (26) 

x 2 x x t) 

- - — If In / + (1 - /) ln(l - /) - / ln(z - 2)]. 

x 

is more accurate than that by Gibbs and DiMarzio. The 
Milchev's expression is recovered from (|25|) setting V = 0, 
which could be easily verified. Thus, our theory is capa- 
ble of capturing the athermal part and the semiflexibility 
part of the configurational entropy at constant pressure 
more accurately than the previous lattice theorie s 2 ' 19 . 
However, we disagree with the conclusions of Ref£ re- 
garding the relevance of the Gibbs-DiMarzio work. Our 
point of view will by disclosed in the next Section. 



V. KAUZMANN PARADOX 

Typical curves of entropy vs. reduced temperature for 
various values of the flex energy, incorporated in the pa- 
rameter Tj = zUf/2, are shown in Fig. [21 As evident 
from the graph, the configurational entropy vanishes at 
a positive temperature T£, and the lower the flex energy 
the smaller this temperature is. It is important to no- 
tice that, for small flex energies, the entropy behaves like 
it has reached a stable level just before the rapid drop 
for T* near TJ . Thus, in an experiment, being unable 
to cool the system down to temperatures below TJ one 
might wrongly conclude that the entropy extrapolates to 
a positive value at T = 0. 

We understand that no state in Nature could by any 
chance have a negative entropy. Therefore, the predic- 
tion of solutions corresponding to negative entropies must 
be an artificial feature of a model and these solutions 
should be replaced by a unique state where the system 
gets trapped under further cooling, i.e. the ideal glass 
state, see£. 

The principal contribution that makes the entropy neg- 
ative is — ( 9x ~ 9 ^ lnr) ; which is one of the terms of the 
thermal correction due to the attractive mer-mer inter- 
action. It remains relatively constant over a wide tem- 
perature range, while other positive contributions reduce 
with temperature decrease. The flexibility term vanishes 
because the polymers becoming stiffer loose their internal 
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= 0.26287 
= 0.264 
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FIG. 2: Configurational entropy dependence on reduced tem- 
perature for different values of Tj calculated taking P* = 
0.25, z = 12, x = 1000. 



degrees of freedom. In the limit of infinitely low temper- 
atures / -»■ and, hence, [/ln(/) + (1 - /) ln(l - /) - 
/ln(z — 2)] — > 0. This suggests that the results of the 
Monte Carlo simulations by Wolfgardt et. aL— cannot 
used to doubt the relevance of the Kauzmann paradox 
due to the absence of the interaction between monomers. 

Another interesting aspect of the glass transition is the 
existence of routes of vitrification alternative to super- 
cooling such as polymerization and compression by pres- 
sure. The first investigation of molecular weight depen- 
dence of the ideal glass transition is contained in the orig- 
inal paper by Gibbs and Di Marzio 2 . Figure [3J presents 
configurational entropy vs. chain length at a constant 
pressure plotted for various temperatures. The curves 
have a common behavior: the entropy is a monotonically 
decreasing function of molecular weight, and the lower 
the temperature the smaller the entropies are. When 
a curve crosses S n = ordinate at ik, further molecu- 
lar wight growth is thermodynamically prohibited, which 
signals the glass transition to occur during polymeriza- 
tion at some x g < xk- Entropy curves descend rapidly 
at initial stages of equilibrium polymerization, then the 
curves becomes almost flat although not bounded by a 
horizontal asymptote. Even a minute temperature dif- 
ference can influence Xk dramatically. This situation is 
clearly reflected in Fig. [3J where the solid and the dashed 
lines correspond to T* = 0.26287 and T* = 0.264, respec- 
tively, and Tj = 1.0 and P* = 0.25 were used. Based on 
this one can introduce an effective temperature above 
which entropy is practically always positive. Calculated 
as accurately as in 5 decimal digits, due to high sensi- 
tivity of entropy, this temperature is T* = 0.26287 for 
T* = 1.0 and P* = 0.25. 

The states with negative entropies are also produced 



FIG. 3: Configurational entropy vs. chain length, P* = 0.25, 
z — 12, Tj — 1.0. This demonstrates the ideal-glass tran- 
sition at high molecular weights and low temperatures that 
corresponds to the system vitrification during polymerization 
reactions conducted at sufficiently low temperatures. 
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FIG. 4: Configurational entropy vs. pressure, x = fOOO, z = 
12, Tf = 1.0. This is the demonstration of the ideal-glass 
transition generated by a combination of low temperatures 
and high pressures. 



at elevated pressures below a certain temperature. The 
shape and the sensitivity of the curves in Fig. [4] are 
similar to the previous case. The entropy is practically 
flat at high pressures and is very sensitive to T*. Taking 
TJ = 1.0 we find that the entropies cross the S n = ordi- 
nate for the temperatures below T* = 0.264. In this case 
the point where the entropy turns negative, P* = P^, 
is the thermodynamic manifestation of the glass transi- 
tion for P* < Pj£. The present theory is in qualitative 
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VII. CONCLUSION 




FIG. 5: Three typical dependencies of the configurational spe- 
cific heat on chain length. For this figure we used TJ = 0.5, 
P* = 0.25, z = 12. 



agreement with the work by Di Marzio et al 
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VI. SPECIFIC HEAT 



One of applications of this theory could be a compari- 
son with experimental results available from literature for 
specific heat dependence on molecular weight. The spe- 
cific heat of linear n-alkanes as a function of the number 
of carbons in the backbone is found to be a monotoni- 
cally decreasing function^, while the specific heat of a 
continuously polymerizing system was measured to pass 
through a maximum^. We suggest that the controversy 
can be resolved in the framework of the present theory 
by studying the configurational specific heat as a function 
of chain length. Figure [5] shows C v vs. x in three cases: 
T* = 0.48, 0.47, 0.46 represented by continuous, dashed, 
and dotted lines, respectively; in all three cases Tj? = 0.5. 
By showing this we want to demonstrate three types of 
typical behavior: for high temperatures (T* — 0.48 and 
above) the curve is always sharply decreasing, then reach- 
ing a constant value; for low temperatures (T* = 0.46 
and below) the curve is monotonically increasing; in the 
intermediate region, the curve initially decreases reaches 
a minimum and then increases. This disagrees with the 
statistical calculation by Wang and Johari^. As was 
shown by Di Marzio and Dowcll 2 t the specific heat has a 
large vibrational part. Hence, the specific heat maximum 
observed by Johari et al£^ during a polymerization reac- 
tion could be due to a freezing of the vibrational degrees 
of freedom as polymers become longer. 



We provide a very careful derivation of the entropy and 
the pressure expressions based on the Bethe-Peierls ap- 
proximation for the model of polymer melt. The predic- 
tion of the system density at a given reduced pressure and 
temperature is shown to be in quantitative agreement the 
experimental data obtained from 1 -. The overall config- 
urational entropy of the polymer fluid consists of three 
fundamentally different types of contributions, i.e., the 
athermal entropy (translational and configurational de- 
grees of freedom) , the semiflexibility correction entropy, 
and the thermal correction entropy that is due to the 
thermal correction term of the partition function. The 
third type, that is absent in the Flory-Huggins theory and 
in the LF theory by Sanchez and Lacombe, contributes 
the main part at low temperatures. We also demonstrate 
that our theory is equivalent to the Ryu-Gujrati theory 
with the correction due to chain semiflexibility . 

The present work suggests that the thermodynamic 
scenario of the glass transition proposed by Kauzmann 
is relevant for the vitrification in the course of a poly- 
merization as well as in the case of the compression un- 
der a sufficiently high pressure. We obtain the negative 
entropy states at certain conditions: low temperatures, 
T C [0..Tk], and some fixed pressure; high pressures, 
P C [Pk---oo), and a fixed temperature; high molecular 
weights, x C [xk-.oo), while both temperature and pres- 
sure are kept constant. These results are in quantitative 
agreement with earlier work by Di Marzio et al&2&. 

Our computation of the configurational specific heat 
as a function of chain length reveled that this function is 
monotonically increasing at low temperatures and mono- 
tonically decreasing at high temperatures. In the transi- 
tion region the function has a minimum. 
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IX. APPENDIX 

In this appendix we demonstrate that the entropy ex- 
pression developed in the paper by Ryu and Gujratii^ 
with added correction due to polymer semiflexibility is 
equivalent to (|25[) . Formulated for multicomponcnt mix- 
tures their theory can easily be adopted to our system 
with minor corrections. One of these is due to discrep- 
ancy in methodology of accounting for the mer-mer in- 
teraction. In their theory this interaction is incorporated 
(as in the Flory-Huggins theory) through the effective in- 
teraction between occupied sites and vacancies with the 
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energy given by 



the following: 



eio = eio - l/2(en + e o), 

where eio, en, eoo are the energies of physical interactions 
for mer-hole, mer-mer and hole-hole pairs, respectively. 
Since only the second is actually present, eoi = 0, en = 
V, e 00 = 0, one gets eio = -V/2, with w = exp(V/2RT) 
being the corresponding Boltzmann weight. Note also 
that q in their symbolic representation is equivalent to z 
in ours - both signifying the coordination number. The 
seminexibility can be incorporated replacing r = q — 1 
in their paper by r/ defined in ([5]) and adding the term 
— N g In Wf to the entropy [see Eq. (f2"8")l below]. Here we 
summarize the results of their theory. The final equations 
resulting from the Bethe lattice iterative technique are 



J 



2/o,i = (^ + Aiyr 1 J/o,iV'"/)/(l + ^Aiyr 1 J/o,i 1 A/): 
1/(1 -6) = 1 + {z/2)xyl iy r X /r f , 
Qi = l + wXiy^y^/rf, 

where Ai = (z/2)x — (x — 1) and the parameters 
Vii 2/0,1, Qi are evaluated from the relations 



2/o,i 



■;(p - 1) + ^/w 2 (p-l) 2 + 4j> /2, Qi = l+wp/y 0A , 



with p = 0Ai/(z/2)(l — 0)x. Below we list the expres- 
sions for the the densities which are numbers of the cor- 
responding species per lattice site: 



b n ,i = 6/x, 6 = (f>„,i(x - 1), 0oo = z(l - 6)/Ql, O i = z(l - 0) - 20 oo , 
011 = 02/2)0-001 -06 , 0i„ = [z/2)6 - 6 , 0o„ = (z/2)(l-0), 

0„=z/2-0fc, 0u=0?„/0u, 0OO=^On/^«' <^01 = 2 0O«01«/0«, 



where nj i and 06 are the densities of polymers and 
bonds, respectively; 0oo, 0oi and 0n are the densities of 
hole-hole, hole-mer and mer-mer nearest-neighbor con- 
tacts; 000' ^oi an d 0ii are the corresponding athermal 
contact densities; 0i u and Ou are the densities of unoc- 
cupied lattice bonds associated with occupied sites and 
vacancies, respectively; U is the density of unoccupied 
lattice bonds. The adimentional pressure, lo = flPvo 
[}3 = l/k B T), given by 

PPvo = -ln(l-0) + (2/2)m(20«/ 2 ) + (z/2)ln(0° o /0oo). 



can be transformed to the form equivalent to (|24| : 



(3Pvq = - ln(l - d) + (z/2) ln[(l - 6)Qi 



(27) 



Finally, we compared the Ryu-Gujrati configurational en- 
tropy corrected to account for seminexibility, 



S n /k B = (0 n ,i ln[(z/2)r;- 2 x/0^i] -(1-9) ln(l - 9) 

+ 0oiln(0oi/0oi) + 0ooln(0o O /0oo) + 0n ln(0° 1 /0 1 

I 



-N g \nw f )/9, 



(28) 



with (|27j) and verified numerically that they are equiva- 
lent. 



1 W. Kauzmann, Chem. Rev. 43, 219 (1948). 

2 J.H. Gibbs and E.A. DiMarzio, J. Chem. Phys. 28, 373 
(1958). 

3 P.D. Gujrati, J. Phys. A: Math. Gen. 13, L437 (1980); 
P.D. Gujrati and M. Goldstein, J. Chem. Phys. 74, 2596 
(1981). 



4 P.D. Gujrati and A. Corsi, Phys. Rev. Lett. 87, 025701 
(2001); A. Corsi and P.D. Gujrati, Phys. Rev. E. 68, 
031502 (2003)^ 

5 P.D. Gujrati, |cond-mat/0309143 

6 M. Wolfgardt, J. Baschnagel, W. Paul, and K. Binder, 
Phys. Rev. E 54, 1535 (1996). 



10 



7 P.J. Flory, Proc. R. Soc. London Ser. A 234, 60 (1956); 
P.J. Flory, Proc. Natl. Acad. USA 79, 4510 (1982). 

8 H.A. Bethe, Proc. Roy. Soc. A 150, 552 (1935); R. Peierls, 
Proc. Roy. Soc. A 154, 207 (1936). 

9 T.S. Chang, Proc. Camb. Phil. Soc. 35, 265 (1939). 

10 P.J. Flory, J. Chem. Phys. 10, 51 (1942); M.I. Huggins, 
Ann. (N.Y.) Acad. Sci. 43, 1 (1942). 

11 K.F. Freed and J. Dudowicz, Adv. Polym. Sci. 183, 63 
(2005). 

12 J. Dudowicz, K.F. Freed and J.F. Douglas, J. Chem. Phys. 
124, 064901 (2006). 

13 A.R. Miller, Proc. Camb. Phil. Soc. 38, 109 (1941); A.R. 
Miller, Proc. Camb. Phil. Soc. 39, 54 (1943). 

14 J.-H. Ryu and P.D. Gujrati, J. Chem. Phys. 107, 3954 
(1997). 

15 J.K. Roberts, Proc. Roy. Soc. A 161, 141 (1937). 

16 I.C. Sanchez and R.H. Lacombe, J. Phys. Chem. 80, 2352 



(1976); R. H. Lacombe and I.C. Sanchez, J. Phys. Chem. 
80, 2568 (1976); 

I.C. Sanchez in Polymer Blends Vol. 1, Ed. D. R. Paul and 

Seymour Newman (Academic, New York, 1978). 

I.C. Sanchez and R. H. Lacombe, Macromolecules 11, 1145 

(1978) . 

C. R. Milchev, Acad. Bulg. Sci. 36, 1415 (1983). 

E.A. Di Marzio, J.H. Gibbs, P. D. Fleming III, and I. C. 
Sanchez, Macromolecules 9, 763 (1976). 

D. Huang, S.L. Simon, and G.B. McKenna, J. Chem. Phys. 
119, 3590 (2003); ibid. 122, 084907 (2005). 

G.P. Johari, E. Tombari, S. Presto, and G. Salvetti, J. 
Chem. Phys. 117, 5086 (2002). 

J. Wang and G.P Johari, J. Chem. Phys. 116, 2310 (2002). 

E. A. DiMarzio and F. Dowell, J. Appl. Phys. 50, 6061 

(1979) . 



